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We study neutrino emission from direct Urea processes in pion condensed quark matter. In 
compact stars witli high baryon density, the emission is dominated by the gapless modes of the pion 
condensation which leads to an enhanced emissivity. While for massless quarks the enhancement 
is not remarkable, the emissivity is significantly larger and the cooling of the condensed matter 
is considerably faster than that in normal quark matter when the mass difference between u- and 
d-quarks is sizable. 
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Q ■ I. INTRODUCTION 

0^ \ Compact stars are born at temperatures T > 10^^ K (about tens of MeV) and rapidly cool down to 10^ K 
'"^ (about a few MeV) within minutes, and then the cooling process is dominated by neutrino emission for about 10^ 
years during which the temperature drops down to 10® K (about a few KeV) P, Q- Since the typical density 
, in a compact star could be several times the normal nuclear density, matter with strong interaction at such high 
rS i' density and low temperature must be highly degenerate and may be in a superfluidity/superconductivity or hybrid 
JL ] state [1, Neutrino cooling processes in these phases may behave differently from that in the normal phase, 

and the observations of cooling rate can provide constraints on the phase structure of dense matter. 
\ Many cooling mechanisms have been investigated in literatures. If a compact star is in the state of nuclear 
'— matter, the cooling is governed by the modified Urea processes and the temperature dependence of the neutrino 
I emissivity is e ~ [tI. Isl. lol. [Toj . When the proton abundance is larg e enough (> 11%), the direct Urea processes 
, can take place and lead to a larger neutrino emissivity e ~ T® [Tll.ll2| When pion or kaon condensation is present, 
■ such a fast cooling can be realized as well [1, [H, [ij, [lB|. In the nuclear superfluid, the quasi-particle spectrum 
QQ of nucleons is gapped and the cooling mechanism is dramatically changed. At low temperatures T ^ with 

' being the critical temperature of nuclear superfluid, the Cooper pairs can rarely be broken by thermal fluctuations, 
\l \ the neutrino emissivity is suppressed by a Boltzman factor e"^/'^ with A being the energy gap [igI. [l7j. When 
\^ ■ the temperature approaches Tc, the effect of pair breaking and recombination results in an emissivity of e 

' for neutrino pair radiation [isj . 
C ■ ■ Compact stars can be in the state of quark matter [l^,!!^], when the baryon density is high enough. The neutrino 
' processes in normal quark matter are similar to that in nuclear matter. For example, the direct Urea processes 
^ , including the /3-decay {d ^ u + e~ + Pe) and the electron capture (m + d + Ve) in normal quark matter 

lead to an emissivity of e ~ T® [2l|, [l^, [2^ , while the modified Urea processes give e ~ [2^ . Quark matter 
. , at asymptotically high density is in the color-flavor locked (CFL) phase [1^ of color superconductivity (CSC). 
. The CSC phases at moderate density are still unclear due to the complicated non-pert urbative effect. Various 
■ ■ ' candidates have been proposed, such as two flavor CSC (2SC) [11], gapless 2SC HJ, gapless CFL (gCFL) [28j . 
crystafline CSC [1^, and spin-1 CSC [13, HU, [111 • Neutrino emission from CSC quark matter is exponentially 
suppressed at low temperature, if the quasi-particle spectra are fully gapped [H, [H, [H, [H, [13, [11|. Similar to 
nuclear matter, the pair breaking and recombination effect results in an emissivity of e ~ at temperature close 
to Tc i39j]. The emissivity behaves as e ~ T^-^ for gCFL [H and e ~ T® for g2SC, see, e.g. Table 1 of Ref. Q. 

At moderate baryon density, the state with spontaneous isospin symmetry breaking is a possible ground state of 
quark matter, if the isospin chemical potential (or equivalently electron chemical potential) is larger than the pion 
mass [40l.l4l|. In this state, the nonzero quark-antiquark condensate (di'f^u) or (ui'f^d) can be considered as pion 
condensate when the coupling between the quark and antiquark is strong enough. When the temperature of the 
system increases, the condensate will melt and there will be, respectively, BCS and BEC phase transitions in weak 
and strong coupling regions. Without making confusion, we call in the following the quark-antiquark condensate 
of light flavors as pion condensate. We emphasize that the pion condensation in the current paper is very different 
from the Goldstone boson one in the CFL phase [H, [4^ where the Goldstone bosons are diquark-anti-diquark 
excitations. In particular the color symmetry is spontaneously broken in the CFL phase while it is kept in our case. 
For example, the kaons K'~' and with quantum numbers ds and us in the CFL phase are actually the lowest 
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excitations and made of {us)(ud) and {ds){ud) respectively, where all diquarks (anti-diquarks) are in anti-triplet 
(triplet) in color and flavor space. In our case the pions are really made of ud, du and uu ~ dd, the same as in 
vacuum. For neutrino emission, the role of pion condensation in quark matter is very different from that in nuclear 
matter. In nuclear matter, the pion degrees of freedom open effectively a new channel for nucleon interaction and 
result in an enhanced neutrino cmissivity [1, In quark matter, however, the pion condensate will suppress, on 
one hand, the neutrino emission of pairing quarks due to the cost of breaking the pairs via thermal excitations, and 
open, on the other hand, more phase space for direct Urea processes of those unpaired quarks. The competition 
between the two effects determines whether the neutrino emissivity is suppressed or enhanced. 

At moderate baryon density, the system is controlled by light quarks and the strange quark excitation is not yet 
important. In this paper, we will calculate the neutrino emissivity and cooling rate through direct Urea processes 
in pion condensed quark matter with only two flavors. Since the typical baryon chemical potential fiB is of the 
order of 1 GeV and the isospin chemical potential fii is about 100 — 250 MeV in compact stars, the mismatch 
between the two Fermi surfaces of the pairing u{u) and d{d) quarks is large enough and the pion superfluid will 
be in gapless state when the condensate A satisfy the condition > A. We will focus on the gapless pion 

condensation. 

The paper is organized as follows. In Section |TT] we derive the general formulation for the emissivity in pion 
condensed quark matter. The calculations of the emissivity and cooling rate for positive energy excitations which 
are dominant at low temperature are presented in Sections IIIII and IIVI The contribution from negative energy 
excitations to the emissivity is discussed in Section [V] The summary is given in the final section. Our units are 
h = ks = c = 1 except particular specification. As a convention, we denote a 4-momentum as K^^ = (fco,k), and 
its 3-momentum magnitude as fc = |k|. 



II. GENERAL FORMULATION FOR NEUTRINO EMISSIVITY 



At a typical temperature for a relatively aging neutron star, the neutrino mean free path is larger than the star 
radius, and there is no accumulation of neutrinos inside the star. In this case, the chemical potential for neutrinos 
can be taken to be zero. We assume /3-equilibrium, so the chemical potentials satisfy fid = fJ-u + /^e for the quark 
chemical potentials fiu and fXd and the electron chemical potential //g. Since the pion condensation favors a high 
isospin or equivalently electron chemical potential, we assume that /Xe is of the same order of fiu and /i^. 




FIG. 1: The one-loop diagram of W-boson polarization tensor. 



It is convenient to use the Fermi effective model of weak interactions to describe the Urea processes where the 
characteristic energy scale is about a few hundred MeV and much less than the mass of W-bosons. The interaction 
Lagrangian of the model is 

Ant = ^J'-Jl (1) 

with the weak current 

J^(a;) = 1^7^(1 -75)e + ?27^(l- 75)d, (2) 



where u^d^e^v denote the fields for u and d quarks, electrons and neutrinos, respectively, and G = Gpcosdc ~ 
1.13488 X 10~^^ MeV~^ is the four-fermion coupling constant with the Cabibbo angle 9c ■ 

The neutrino emissivity, defined as the total energy per unit time and per unit volume carried away by neutrinos 
and anti- neutrinos, can be written as 



/(^3p r ^3p 

xL^''{P,,Pe)lmnf^{Ee - - E,,p, - p,), (3a) 

xL^^iP,, Pe)lmn^,{Ee - + E,, + p,), (3b) 

where Pi — (Ei,pi) for i — w, d, e, are on-shell 4-momenta for quarks and leptons. The energies are Ei = 
Ep. = \/ Pi + rn^ with the masses for i and with = and mg w 0. Here nsix) = (e^/-^ — 1)^^ and 



3 



np(x) = [e^l^ + 1) ^ are the Bosc-Einstein and Fcrmi-Dirac distribution functions, respectively. L^'^{Pe,P,j) 
denotes the leptonic tensor and nj^^(Q) the retarded polarization tensor for W-bosons. The leptonic tensor reads 



L^''{P.,Pe) - Tr[7^(l - 7.5)7 -^eT'll- 75)7-^.] 
As illustrated in Fig. [l]the M^-boson polarization tensor can be written as 



(4) 



(5) 



where prfo = PuQ+Iq, Vd = Pu+q, and = 7'^(1~75)''"± with the ladder operators in flavor space t± — (Ti±ir2)/2. 
The bosonic and fcrmionic Matsubara frequencies are given by qo = i2lT:T and p^o — i{2n+ 1)ttT (/, n are integers) 
respectively. The retarded polarization tensor 11^ is obtained from the above after taking the analytic extension for 
the Matsubara frequency qo = i2lT:T qq + iO^ , where the second qo is real. We can then extract the imaginary 
part Imll^ from 11^. The full propagator for quarks 



S{K) = 

in flavor space can be derived from its inverse [Zlj . 



Suu{K) Sud{K) 
Sdu{K) Sdd{K) 



(6) 



(7) 



where m is the average effective quark mass, and ~ (dij^u) and ~ {ui"/^d) are the charged pion condensates. 
Without loss of generality, we assume real condensates and let A+ — A^ = A. In this case, the four matrix 
elements of S are explicitly expressed as 



Suu — ^ ^ Bj, (fc) 



a,r— lb 



Sud 
Sdu 



ko + fiB/S + rE^ 
Ar7o 



a.,r— ± 

-t ^ r^Bl{k)B-_{k) 

a,r— lb 



A^75 



a,r— lb 



^ rJ B%{k)Bt{k) 



ko + fiB/S + rE^' 

Afc°75 
ko + fin/i + rE^'' 



(8) 
(9) 
(10) 
(11) 



where the summation over a (and b in the following) is for positive and negative energy excitations for quarks, and 
the that over r (and s in the following) is for quarks and anti-quarks. The Bogoliubov coefficients B^ik), energy 
projectors AJJ and energies E^ are defined by 



Et 



1 — ar 



Ek + aSf^L 



El 



1 + a 



7o(7-k-|-m) 



Ek 

^{Ek + aSfi)^ + A^. 



(12) 



We use Sfi = {^d — tJ-u)/2 = — /i//2 = /ie/2 and fi = to replace and fie, where fJ.i.B,e are chemical 

potentials for the isospin, baryons and electrons respectively. From the poles of the quark propagators, the quasi- 
quark energy is then given by E^ + rfi. The chemical potentials for u and d quarks can be expressed in terms of 
Sfi and fi ds. flu = H — 5fi and = ^ + (S/i. 

Substituting the quark propagator ^ into the W-boson polarization tensor ([S]) and performing the trace in 
color and flavor space, we obtain 

n^'^(go,q) ^ NeTj2j ^Tr[7^(l-75)S,„(p„o,Pn)7"(l-75)^dd(pdO,Pd)] 



E 



a.b.r,s—:L 



d^Pu rr\a,r. r. ^ r.a , ^ r.~b , F{~r E^^ - fl) - U p {,~ sE ~ fl) 



{2'KY'iEuE, 



■H^S{Pu,Pd)B-r{Pu)BZl{pd)- 



qo - rEi + sE- 



,(13) 
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with the quark tensor H^^{Pu, Pd) given by 

H^biPu,Pd) = iEuEdTT[^\l - 75)A;^ 707^(1 - 75)A^,7o] 

= 8{PauPbd + PLPL - Pan ' PmQ^'' + »e V/J^au^'^d] > (14) 

where Pau = (E'uiaPu) with a = +, — and Pm = {Ed,bpd) with 6 = +, — for positive and negative energy 
excitations respectively. We have also used the number of colors Nc which is finally set to 3 in our calculation. 
Using Imn^°'(qo, q) = linIl'^'^{qo + iO+, q), the imaginary part of the retarded polarization tensor of W-bosons can 
be read out directly, 

Imn^'^(go,q) = ttN, Yl J (27rffigl£:, ' iPu)BZ':{pd)ns\~qo) 

xnpi-rE^ - fi)nFisE-'' + fi)6{qo - rE'i^ + sE^'') 

a,6,r,s— lb 

xnpirE'^ + y)nF{'sE-'' - ^i)S{qo - rE^ + sE-"), (15) 

where the first equality is for the neutrino emissivity with qQ = E^ — He ~ Ey and the second for the anti-neutrino 
one with go = ^'e ^ Me + Ei,. Here we have used the identity np{x) — npijj) = nF{—x)np{y)/nB{—x + y) = 
-nF{x)nF{-v)/nB{x - y). 

Substituting the imaginary polarization (jlSp into the neutrino emissivity ([3]), we arrive at 

a.h,r,s—± 

x53(p^ - + p„ - Vd)B^{v^^B-l{Pd)nF{Ee - fie)nFi-rE^ - ti)nF{sE-^ + n)\Mab? , 

/ i2n)^2Ej {2^Y2eJ {2nY2Eu J {2nY2Ed > ^ - " d ) 

x53(p^ + p, + p„ - Y*d)B^r{Pu)BZl{pd)nF{t^e ^ E^WirE^ + ^i)nF{-sE^'' - Ai)|M,b|', (16) 

where the shorthand notation |A/af,p is the spin- averaged matrix clement for the /3-dccay or the electron capture 
processes [2lj |. 

|M,fcp = ^LxAPe,P.)H^biPu, Pd) = GiG^iPe ■ Pau){P. ' Pbd) ■ (17) 

Here we have suppressed energy projection indices a, 6 in quark momenta as before. When a = + and 6 = + 
for u and d quarks respectively, |M++p is just the matrix element of the Urea processes in vacuum. For other 
combinations of a and h, it represents the matrix elements with at least one negative energy excitation for u, d 
quasi-quarks. 



III. NEUTRINO EMISSION WITH POSITIVE ENERGY EXCITATIONS 



It is well known that the negative energy excitations are normally neglected at low temperature due to their 
strong suppression relative to the positive excitations, see, e.g. Eqs. (21-24) for fully gapped phases in Ref. [s^ . 
This can be understood from the analysis of the parts involving Bogoliubov coefficients and quark number distri- 
butions in Eq. ([16]), 

B'^{p^)Bz'i{Pd)nF{TrEl T pi)nF{±sE-'' ± fi), (18) 

where the upper and lower signs correspond to neutrinos and anti-neutrinos respectively. Let us consider the 
part related to the negative energy excitations for u quarks, B~ {pu)nF{—rE~ — /i) and B~ {pu)nF{rE^ + h). 
Note that the dominant contribution to the emissivity comes from the phase space near gapless nodes which 
correspond to the vanishing arguments of quark distribution functions by selecting the branch with r = — , see 
Appendix [X] The gapless node for u quarks in this case is located at E'" = ■\/ /.j^ — + then we obtain 

BZ{pu) = ^1 — i/ j-j? — A2//i^ /2 which approaches zero for large fi or small A. The same conclusion is also true 
for d quarks. The above arguments are valid only at low temperatures T <^ fi^ d where the dominant contribution 
of phase space integrals comes from the region near gapless nodes. 

In this section we calculate the neutrino emissivity by dropping out the negative energy excitations. Wc will 
briefly discuss the contribution from the negative energy excitations in Section |Vl 
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When quarks are treated as free and massless constituents in normal quark matter, phase space for direct 
Urea processes vanishes, since quarks and electrons are collinear in momenta on their Fermi surfaces plj . To 
get non-vanishing phase space which requires that the Fermi momenta of u and d quarks and electrons should 
satisfy a triangular relation, one must take into account nonzero quark masses and/or modified quark dispersion 
relations due to interactions and/or pairings between (anti-) quarks. We will consider these effects together. When 
the baryon chemical potential is high enough, the dispersion relations of quasi-quarks are gapless, the neutrino 
emissivity is dominated by the quasi-quarks around the gapless momenta. For u and d quarks, the gapless momenta 
are = {I — (i?,")^ — and = (1 — k)-^(£'^)2 — with the gapless energies E'^ = _ ^2 „ 
and = y/jj^^^A^ + 5^. The dispersion relations of quasi-quarks are schematically shown in Fig. O We have 
included the interaction among quarks by introducing a phenomcnological constant k which reflects the reduction 
of the Fermi momenta and can be derived in Landau's Fermi liquid theory or in flcld theory. In perturbative 
QCD we have k = 2Q;s/(37r) with ag being the strong coupling constant 4J], and in the NJL model it can be 
expressed as k = 4(7/i^/(37r^) with the four-fcrmion coupling constant g (38|. In the normal phase with vanishing 
pion condensate, E^ ^ and ^ arc just the Fermi energies and momenta, i.e. E'^ ^ — ^u,d and Pu d — '^)f^u,d- 



2Q0 




IPl (MeV) 



FIG. 2: The schematic dispersion relations for gapped and gapless quasi-quarks. 



At the gapless momenta, the angles between u-quark and electron momenta and between u- and d-quark 
momenta can be written as 

Substituting them into Eq. (|17p . one obtains approximately 



Af++|2 « 6AG^EeE.,E,Ea (l-^ cos^ ) ( 1 - cos(?d, ) , (20) 



EO «V V EO 

where 9dv is the angle between neutrino and d-quark momenta. In the brackets the u- and d-quark energies and 
momenta and the angle between u-quark and electron momenta are replaced by the corresponding values at the 
gapless nodes. Taking into account the relation /.td = ^lu + /-te, the (5- functions in Eq. (jl6p can be simplified as 
^jLe/ {PuPd)S{cos0ud — cos0^^) to the leading order of p^/^^, since the higher order contribution is small in the 
kinetic range of the processes. Using the above approximation, the neutrino and anti-neutrino emissivity can be 
expressed in the following form 

« 16^iVeGVe-B°i^S 1- f^cosC M I] / dE.El dEu dEd 
xB^ {pu)BZs{Pd)nF{Ey + rE+ - sE^)nF{-rE+ - ^i)nF{sE^ + fi), 

/ \ roc POO POO 

« 16^iVcGVei^°i^S 1- l^cosC M I] / dE.El dEu dEd 

xB+{Pu)BZMnF{E, - rE+ + sE;i)nF{rEt + ii)nF{~sE^ - /i), (21) 
where we have used dp^dpd = dEudEdEuEd/{puPd)i and the angular integration J has been evaluated as 

dO,,, f dflu f dQd Pd „ „ „n , 1 



^ ^ y (# y (2# y (2^(1 - ^cos...).(cosa.. - cosC) = (22) 

To further simplify the emissivity, it is convenient to introduce the dimcnsionless variables 

Eu + S^i Ed~ Sn ^ E^ A 

,y= j; , V ^ —, w = (23) 
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and relax the integration ranges of x^y to (—00, cx)), since the integrations are dominated by the gapless nodes. 
Dropping out odd terms of a; and y from the integrations, and noting that the summation over r and s is symmetric, 
and ep in ([2T|l are approximately the same. Therefore, we have 



AGVe^;X°(l"f COS( 



T^F{u,w), 



(24) 



with 



F{u,w) = ^ Frs{u,w), 



r,s— lb 



1 



1 



Frsiu.w) = / dxdydvv^ , ., , , 7^^^^ ■ 



(25) 



While the two integrals in (|2ip for neutrinos and anti-neutrinos are different in general case, see Appendix [B] they 
approach to the same value at low temperature, because the gapless momenta scaled by temperature are far from 
zero providing high effective Fermi surfaces. In this case, the integrals perfectly center at the gapless momenta 
and other contributions quickly damp out. However this is not the case at Tc- The difference and coincidence in 
the neutrino and anti-neutrino emissivities also occur in fully gapped superfluid/supeconducting states. 

Eq. (|24p is our main result, which presents the neutrino emissivity e as a function of temperature T, chemical 
potentials hb and /Xe, the effective quark mass mu,d, and the pion condensate A. 

To clearly see the effects of the pion condensate, we consider two limit cases. One is the high temperature 
limit with T close to Tc where the condensate is small, A <C fj,u, Hd, fJ'ejT, and its temperature dependence 
can be modeled as A = Ao^l — {T/TcY with Aq being the pion condensate at T = 0. The other is the low 
temperature limit with T <^ Tc. Since the chiral symmetry is largely restored in the phase of pion condensation, 
the effective quark mass is small compared with the quark chemical potential (actually, 7Ti//i cx y^), we can assume 

TO < flu,fJ.d,l^e- 

In the hmit of high temperature, we can expand ^, p° ^ and cos0"g in terms of k, to^ ^//i^ and A^//x^, see 
Appendix ICl By keeping only the leading order and considering the limit Limyj^oF {u,w) = 45 77r^/5040 [i^, we 
obtain the emissivity near Tc, 



1 



A2 



1371 



Hu^idJ 2520 



TtG fJ.efJ'ufJ'd 



fJ-u 



(26) 



Comparing it with the result in normal quark matter with massless quarks [2l|, 

1371 



2520 



^gVmum4 1-1^ cos Cjr' 



E" 



1371 2 , \ rr.6 

2520"^ MeMuM. (^1 + - Ur 



(27) 



we conclude that, near Tc a small pion condensate reduces the emissivity slightly but the correction due to the 
quark mass splitting may enhance it significantly. 

In the low temperature limit, at = (w = 0) we can recover the fully gapped phase and the emissivity has an 
exponential suppression factor e~^l'^ (see, e.g. |35|). For non- vanishing [Ib or m, such a suppression is absent due 
to the gapless modes appeared in quark number distribution functions. In the case with m > w 3> 1, we derived 
in Appendix [Al the asymptotic behavior of Frs{u,w) which can be summarized as 



F++ 


- we"" 


F+- 




F-+ 






- 0(1). 



in' 



(28) 



It is easy to see that the term F^^{u,w) dominates F{u,w) due to the gapless modes in both u and d quark 
number distribution functions. Therefore, at low temperature the emissivity is of the same order as that in normal 
quark matter. 

To do numerical calculation, we choose the parameters /is = 900 MeV, ni = 150 MeV, k = 2/{3tt) and 
Aq = 100 MeV. These are typical values for the possible pion condensed quark matter in compact stars. The 
numerical result for the ratio e/eo as functions of T/Tc is illustrated in Fig. [3] where the critical temperature is 
given by Tc = 0.57Ao according to the ECS theory. When the chiral symmetry is fully restored, the quarks are 
massless, and the difference between pion condensed and normal quark matter is small at low temperature and 
becomes considerable at intermediate temperature. At high temperature T Tc, the emissivity in pion condensed 
matter is a little lower than that in normal matter, just as we expected above. However, in general case with 
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Ao=100MeV 



mu=75MeV 
md=125MeV 



mu=md=OMeV 



0.2 0.4 0.6 0.8 1 

T/Tc 

FIG. 3: The scaled neutrino emissivity e/eo as functions of the scaled temperature T /Tc for massless and massive quarks. 

nonzero baryon and isospin chemical potentials, the effective quark masses m„ and for u- and d-quarks are 
not the same, there should exist a quark mass splitting induced by the isospin symmetry breaking. In the massive 
case with rriu = 75 MeV and md = 125 MeV, the enhancement factor for the emissivity becomes 1.2 in the whole 
temperature range. A similar behavior was found in a two-flavor gapless color superconductor [36l |. 

IV. COOLING RATE 



In this section, we compute the cooling rate due to neutrino emission. To this end, we integrate the energy 
equation e{T) = —cv(T)dT/dt and obtain 



t-tn 



.,cv{r) 



dT , , 

To <T') 



-2^= E / -^{nF{El+r^,)\nnp{Et 



where Tg is the initial temperature at time tp- From the entropy density 

a,r=± 

+ [1 - nF{El + r^l)] ln[l - np{El + r/x)]|, 

the specific heat cy is given by 

.ds 



rfi) 



cv{T) 



T- 



dT 



E P'E'a 



(29) 



(30) 



(31) 



where we have used the relation A(T) = Aq^I — {T/TcY and made the approximation that the integral is 
dominated by the gapless nodes, and the functions K{u,w) and G{u,w) are defined by 



K(u, w) 
G(u, w) 



/>oo 

E / ^^(^ 



+ u)2 + ru) 



r=± 



E 



dx 



(32) 



Note that the contribution from the excitation branch with a = + is suppressed by a factor e~^^'^ at low temper- 
ature, it is negligible except near the critical temperature Tc- At Tc with w = 0, we recover K(u,0) ~ 7r^/3 in 
the normal phase. The first term of G'(w, 0) is approximately unit and leads to the standard jump of specific heat 
at the critical temperature for a symmetric superfluid without Fermi surface mismatch between the two species. 
The second term of G{u, 0) reflects the gapless nature of the asymmetric pairing. If one turns off the interaction 
between quarks {k = 0), the X-tcrm reproduces the low temperature specific heat of two- flavor quark matter 
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Cyo = iT with 7 = E^p^ + E^p^ with the u- and d-quark Fermi energies and Fermi momenta, while the G-term 
gives the jump of the specific heat at the normal-superfluid transition point. The numerical result of cv{T) is 
shown in Fig. ID The quark mass dependence of cy is weak and can be neglected. 
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FIG. 4: The scaled specific heat cv/{'yTc) as functions of scaled temperature T/Tc. The straight dashed line is the result 
in normal quark matter, and the solid line is the result in the pion condensation. There is a jump at the phase transition 
point. 



Substituting the neutrino emissivity ([24)) and the specific heat ([3T|) into Eq. (|29|) . we obtain 



t-tn 



dT plEl+p^^E^^ K{u,w) + ^lT-^G{u,w) 



F{u, w) 



(33) 



where ^ and E'J] ^ depend on T through A(T). Fig. [5] shows the time evolution of the pion condensed quark 
matter with initial temperature Tq = 1 McV at tg = 1 yr. For comparison we also show the cooling evolution 
in normal quark matter. In the case with mass difference between u- and d-quarks, the cooling of the condensed 
phase is faster. This is due to the larger neutrino emissivity and smaller specific heat in the pion condensed quark 
matter. For massless quarks, the two cooling curves for condensed and normal quark matter are almost the same 
due to the fact that the emissivity difference between the two cases is small at low temperature, see Fig. |31 




> 60 
lU 



Log(t/yr) 



mu=md=0 



\. normal 

Ao=100MeV~^ 



Log(t/yr) 



FIG. 5: The cooling curves in pion condensed (solid lines) and normal (dashed lines) quark matter for massive (left panel) 
and massless (right panel) quarks. 



V. CONTRIBUTIONS FROM NEGATIVE ENERGY EXCITATIONS 



As we have argued above that the positive energy excitations dominate the neutrino emissivity at /ig 3> A and 
at low temperature T <^ fJ.u^d- In principle, the contribution from negative energy excitations must be considered 
if we follow the exact calculation p6|) at high temperature. In order to better understand the role of negative 
energy excitations, we discuss the emissivity at T — * Tc- In this case, the j, k, b, c related part of the integrands in 
Eq. pB]) reads 

= B!^{p.u)Bz':ipd)nF{E, + rE^ - sE-'')nF{~rE!^ - fi)nF{sE-'' + fi)\Mab\^ 

~ e[-raSgn{Eu + aSfi)]e[-sbSgn{Ed - bSfi)]nF{E^ + r\E.u + aSfi\ - s\Ed - b5fi\) 

xnF{~r\Eu + aSfi\ - ^,)nF{s\Ed - bSfi\ + fi)\Mab\'^ 
= HFiE^ - aEu + bEd - pLe)nF{aEu - Hu)nF{-bEd + Aid)|A/ab|^, 
h — nFiE^ + aEu - bEd + fJ.e)nF{-aEu + fJ,u)nFibEd ~ fJ,d)\Mab\'^, (34) 
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where we have used the approximation B^{pu) ~ 6 [—raSgn{Eu + aSu)] at T ^ Tc which means that only the 
terms with rSgn(£'„ + aSfx) = —a are relevant in the calculation of neutrino emission, the step and sign functions 
are defined as Q{x) = 1 for x > and for a; < 0, and Sgn(x) = 1 for x > and —1 for a; < 0. 

The above expressions with a = + and b = + reproduces the conventional result for normal quark matter, 

h - nF{E^ + Eu- Ed + fie)nF{-Eu + fi.u)nF{Ed- Hd)\M++\^- (35) 

At low temperature T <C this contribution dominates the neutrino emission, since the integrals in p6p have 
sharp peaks at E^ ^ ju„ and E^ ~ fi^ and the contribution from negative energy excitations is suppressed at least 
by a factor of e"'^" ''/^. However, at temperature T ^ l^u.dj such a suppression disappears, and those processes 
which are forbidden by energy conservations at low temperature can happen. For example, for ab = — h the 
processes with energy conservations E^ + Eu + Ed = He and Eu + Ed = E^ + /Xe for neutrinos and ant-neutrinos are 
forbidden at low temperature, but allowed at T ~ ^„_(j since the quark energies could be as large as the chemical 
potentials. 

VI. SUMMARY 

Quark matter may exist in the core of a compact star. Pion condensed phase is one possible ground state of the 
quark matter due to the non-zero isospin density as a result of the (3 equilibrium between u and d quarks. For the 
neutrino emission, the role of the pion condensation in quark matter is different from that in nuclear matter. In 
nuclear matter, the presence of pion degrees of freedom eff'ectively opens a new channel for the nucleon interaction 
and results in an enhancement of the neutrino emissivity. In this paper we considered the gapless pion excitations 
which can be realized in compact stars. We calculated the neutrino emissivity and cooling rate through direct 
Urea processes in pion condensed quark matter with two flavors. We derived general expressions for neutrino 
and anti- neutrino emissivity. At temperature much lower than quark chemical potentials, we demonstrated that 
the positive energy excitation dominates the emissivity if the baryon chemical potential is much larger than the 
magnitude of the condensate. We also discussed the role of negative energy excitations at temperature comparable 
to the chemical potential. 

We showed the numerical results for the neutrino emissivity and cooling rate at low temperature. The total 
neutrino emissivity is of the same order as that in the normal phase, which is a result of gapless excitations. While 
the results in the two phases are quite close to each other when quarks are massless, the mass difference between u- 
and d-quarks due to isospin symmetry breaking leads to a remarkable enhancement of the emissivity. For typical 
values of the initial temperature, the chemical potential and the pion condensate corresponding to the compact 
star interior, the cooling of the condensed phase is found to be faster than that of the normal one, when the quark 
mass splitting is sizable. 
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the startup grant from University of Science and Technology of China in association with the 'Bai Ren' project of 
Chinese Academy of Sciences. 
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APPENDIX A: ASYMPTOTIC BEHAVIOR OF F{u, w) AT LOW TEMPERATURE 



In the limit of u > w 3> 1 at low temperature, the four components of F[u, w) can be simplifies as 

-1 



dvdxdyv 



dvdxdyv 



we 

coo 



F_| {u,w) ~ / dvdxdyv 



dvdxdyv^ e-'^^^+^~'' eV^^^ + e" 







dxdydvv 



1 

4 

1 

4 



1 



F--{u, w) 



where we have used 



dx / dye-V^^~" 

U;l/2(u2-«;2)9/2g-«-^«^ 

/•oo 

/ dv 



\/x^ + + y/y^ + 



dx 



\^u^ — w^ J 



o(i), 



Lim^^o- 



e(-a;). 



(Al) 



(A2) 



The dominant contribution to the integral over x in comes from the regions near the upper limit x 

\/ — vp- — 0^ and lower limit x = \J v? ^ vfl + 0+. 



APPENDIX B: EXACT CALCULATION OF EMISSIVITY FOR NEUTRINOS AND 

ANTI-NEUTRINOS 

In this appendix, we calculate exactly the neutrino and anti- neutrino emissivity listed in (|21|) . in order to see 
how good is the approximation we made in deriving (|24p . From (j2ip . the pre- factors in and Cp are exactly the 
same, and the difference is from the summations and integrals. We denote the summation and integration in ty 
by Fv and that in ep by Fp. It is easy to see the relation F„{—Sfj,) = Fi^{Sfj,). In fact, it is true even for the original 
Ej/ and ep in ([T6|) . For fi ^ 6^, the two integrals are dominated by the phase space near gapless momenta at 
low temperature, the dependence of F^ and Fp is negligible, and we have therefore Fi, = Fy. However, if the 
temperature is not very low or 5^ approaches to ^, the gapless regions might not dominate the integrals, and the 
exphcit dependence on 5^jl must enter and F^. In Fig. [SI the exact F^, and i^p scaled by Fq = -F^(Ao = 0, T = 0) 
arc presented as functions of T jT^. While at low temperature F^, and Fy coincide very well, indicating that the 
approximation we used is good, the difference between F^, and F^ increases as the temperature grows and reaches 
10% at T -> Tc. 
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FIG. 6: The scaled -F^ and Fj? as functions of scaled temperature T/T^. 



APPENDIX C: PARAMETER EXPANSION NEAR Tc 



At T —> Tc in the pion condensation, to the leading order of k, w? / fj? , A^//i^, the gapless energies and momenta 
of quarks and the corresponding angle between w-quark and electron momenta can be written as 



Pd 



A2 



A2 



fJ-u 



I- K 
I- K 



COSf 



which lead to 



E^E^, \ 1 



E?. 



COS 61° 



f^uPd [ 1 
X 1 



1 - K- 
fJ-ul^ld 



A 



2 1 



2/4 2^„/i 
m?, A^ 



2/.tdM 



'2l-iHd J 



1 - K 



flu 2/i2 2/i„/ie 



2/iu/^e 



(CI) 



(C2) 



Substituting this relation into we obtain (|26p . 
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